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Abstract 


The static electromagnetic fields are studied here based on the standard spaces of the physical presentation, and the 
modal equations of static electromagnetic fields for anisotropic media are deduced. By introducing a set of new poten- 
tial functions of order 2, several novel theoretical results were obtained: The classical potential functions of order 1 
can be expressed by the new potential functions of order 2, the electric or magnetic potentials are scalar for isotropic 
media, and vector for anisotropic media. The amplitude and direction of the vector potentials are related to the anisot- 
ropic subspaces. Based on these results, we discuss the laws of static electromagnetic fields for anisotropic media. 
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1. Introduction 


By the Maxwell’s electromagnetic field equations, we 
know that the electric and magnetic field are independent 
each other under the condition of static fields. The clas- 
sical electromagnetic field theory also believes that the 
static electric field can be described by a scalar potential 
function, and the magnetic field by a vector one. Fur- 
thermore, for the passive region, the magnetic field can 
also be described by a scalar potential function [1,2]. 
But it should be pointed out that these results can only be 
obtained in the condition of isotropy, and are also only 
suit for the isotropic media. However, with the develop- 
ment of material science, more and more anisotropic 
dielectric or magnetic materials are applied to various 
fields, such as electron devices, communications and 
sensors, even for the traditional geological structure, we 
also can see the electrically anisotropic media or mag- 
netically anisotropic media. It is found by recent re- 
search works that the limitations of classical static elec- 
tromagnetic field theory have become obvious for these 
anisotropic media. For example, the above results for 
isotropic media don’t exist for anisotropic media, even 
we don’t know the definite form of the electric field po- 
tential function or magnetic field potential function, 
which make a great difficulty in solving the problem of 
anisotropic static electric or magnetic fields [3-5]. Unlike 
the classical static electromagnetic field theory, which 
studies the Maxwell’s equations under the geometric 
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representation, in this paper, the Maxwell’s equations are 
restudied under the physical representation. As the result 
of this, the modal equations of static electric or magnetic 
fields are deduced, which give the novel expressions for 
the potential functions of static electric or magnetic 
fields for anisotropic media, and bring to light the intrin- 
sic laws of static electromagnetic field. 


2. Standard Spaces of Electromagnetic 
Media 


In anisotropic electromagnetic media, the dielectric per- 
mittivity and magnetic permeability are tensors instead 
of scalars. The constitutive relations are expressed as 
follows 


D=exE (1) 
B=uxH (2) 


where the dielectric permittivity matrix € and the 
magnetic permeability matrix x are usually symmetric 


ones, and the elements of the matrixes have a close rela- 
tionship with the selection of reference coordinate. Sup- 
pose that if the reference coordinates is selected along 
principal axis of electrically or magnetically anisotropic 
media, the elements at non-diagonal of these matrixes 
turn to be zero. Therefore, Equations (1) and (2) are called 
the constitutive equations of electromagnetic media un- 
der the geometric presentation. Now we intend to get rid 
of effects of geometric coordinate on the constitutive 
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equations, and establish a set of coordinate-independent 
constitutive equations of electromagnetic media under 
physical presentation. For this purpose, we solve the fol- 
lowing problems of eigen-value of matrixes 


(e-Al) 7 =0 (3) 


(u-yI)p=0 (4) 
where A, (i =1, 2,3) and g; (i =1, 2,3) are respectively 
eigen dielectric permittivity and eigen magnetic perme- 
ability, which are constants of coordinate-independent. 
j,(i=1,2,3) and g(i=1,2,3) are respectively eigen 
electric vector and eigen magnetic vector, which show 
the electrically principal direction and magnetically prin- 
cipal direction of anisotropic media, and are all coordi- 
nate-dependent. We call these vectors as standard spaces. 
Thus, the matrix of dielectric permittivity and magnetic 
permeability can be spectrally decomposed as follows 


6=@0AQ' (5) 
u= YI" (6) 
where A=diag|/,,4,,4,] and TI= diag|7,72:7] 


are the matrix of eigen dielectric permittivity and eigen 
magnetic permeability, respectively. @ = { Jie js} and 


Y ={0, p,p} are respectively the modal matrix of 


electric media and magnetic media, which are both or- 
thogonal and positive definite ones, and satisfy 
OOa ly Pe aT:, 

Projecting the electromagnetic physical qualities of the 
geometric presentation, such as the electric field intensity 
vector E , magnetic field intensity vector H , magnetic 
flux density vector B and electric displacement vec- 
tor D , into the standard spaces of the physical presenta- 


tion, we get 
D'=j7xD_ i=1,2,3 (7) 
E; =j'xE i=1,2,3 (8) 
B'=ọ7xB_ i=1,2,3 (9) 
Bap i=1,2,3 (10) 


These are the electromagnetic physical qualities under 
the physical presentation. 

Substituting Equations (7)-(10) into Equations (1) and 
(2) respectively, and using Equations (5) and (6) yield, 
we have 


{D} =[A]{z} (11) 
{2} =[]{"} (12) 


or 
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D? =4E  i=1,2,3 (13) 
i=1,2,3 (14) 


The above equations are just the modal constitutive 
equations in the form of scalar. 


B=y,H, 


3. Matrix form of Static Electromagnetic 
Field Equation 


The classical static Maxwell’s equations in passive re- 
gion can be written as 


V x E=0, V- D=0 (15) 
Vx H=0, V - B=0 (16) 


where V is a Hamilton operator. It is seen from the 
above equations that the electric field and magnetic field 
are not only independent, but also the same in the form 
of equation. So, it is undistinguishable to study the prob- 
lems of electric field or magnetic field under the static 
condition. For this purpose, we consider here only the 
problem of electric field. 

From Equation (15), we can see that one is a vector 
equation, another is scalar one. It is well known that the 
vector equation can be written as the matrix one, but the 
scalar equation can not. By the first one of Equation (15), 
we have 


[o]{z}=0 (17) 
where 
0 -d, 0, 
[ele | 3. 0 -ə, (18) 
-d, ð, 0 


It is an operator matrix of order 1. 
In order get the matrix expression of static electromag- 
netic equations, both of Equation (15) should be re- 
formed in a suitable form. 

For dynamic electromagnetic fields, a matrix equation 
of electromagnetic waves be dedued by author [6] 


[AE = -v [a]l]; (19) 
where 
-(3 +3) ôy ô 


[4]=[2][@]=]9,. -(2.+0.) 2; CO 
On 


P 
ô, (3n +3) 

It is an operator matrix of order 2. For static electro- 
mag-netic fields, we have 


[A] {£} =0 (21) 


Now, rewriting the second one of Equation (15) in the 
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index form of tensor 
D, =0 (22) 
Differentiating the above equation with index j, it 
become a vectorial one 
D, = 9 (23) 
Rewritting it in the matrix form, we get 
[v]{D} =0 (24) 
where 
Oi On Ox 
[VS ôn ĉn ĉn 
Oe ôn ôy (25) 
In this paper, [A] and [v] is defined as the matrix 
of electic intensity and electic displacement operators 


respectively 


4. Eigen Equations of Static Electric Fields 


Now, we transform the matrix equations of static electric 
field into modal ones. 

Substituting Equation (7) into Equation (21), and mul- 
tiplying it with the transpose of modal matrix in left, we 
have 


[e] [A][@]{z"} =0 (26) 


It be proved [6] that the matrix of electric intensity 
operator can also be spectrally decomposed, that is 


[oy [alle] [a] (27) 
Thus, Equation (26) can be uncoulped and become 
[a ]{z"} =0 (28) 
or 
A E'=0 i=l,2,3 (29) 
in which 
ay =a a) i=1,23 (30) 
where 
(a) = {47 PHA) G1) 


In same way, substituting Equation (7) into Equation 
(24), and multiplying it with the transpose of modal ma- 
trix in left, we have 


[e] [V][@]{D"} =0 (32) 
let 


[oy ely] (33) 
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and substituting Equation (11) into Equation (32),we 
have 


[v |[A]{z"}=0 (34) 
Comparing Equation.(34) with Eq.(28) ,we get 
Lv |[A)=[4] (35) 


It is seen that [v] is also a diagonal matrix. We call 


it as eigen matrix of electric displacement operator. Thus, 
we have 


VD. =O. i=l1,2,3 (36) 


So, Equations (29) and (36) constitute of the eigen 
equations of static electric field. Different from the clas- 
sical ones, they show the simplicity and symmetry of 
static electromagnetic law. 


5. General Solution of Eigen Equations of 
Static Electric Fields 


Let 


Ef =V;%; i=l1,2,3 (37) 


D' =A, i=1,2,3 (38) 


where { x} is an unknown row vector, which is new 


electric potential function of order 2. 

Substituting Equations (37) and (38) into Equations 
(29) and (36) respectively, a unified equation are ob- 
tained as follows 


ox =0_ 721,23 (39) 
where, o,=A,;-V,; (i=1,2,3) is ith modal operator of 
electric field, and a differential operator of order 4. Def- 
erent from the Laplce’s equation for the classical electric 
potential function of order 1, the new electric potential 
function of order 2 can be solved by the modal differen- 
tial equation of higher order, and the classical electric 
potential function of order 1 can be expressed by the new 
electric potential function of order 2. Once the modal 
potential functions are solved from Equation (41), the 
electric intensity and electric displacement can be ob- 
tained by the following conversion 


{E}={a}Vixn, +{B}Vi 1 +i} V2 (40) 

{D} = {PI AV iL + {to} A.V ox + {Os} AsV 325 (41) 
In order to get the classical electric potential function 
of order 1, we rewritting Equation (40) by using Equa- 
tions (30) and (35) 


(EEE {er} x= EEE x) an 
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let 
yaa {ay x (43) 


It is just the electric potential function of order 1 for 
anisotropic media. Thus the electric intensity and electric 
displacement can be expressed by the electric potential 
function of order 1 as follows 


{E} =-{0i vi {23} v2 — {Os} Ms (44) 
{D} =-{d Aw, -{a}Aw.-{GhAw, (45) 


6. The Modal Boundary Condition of Static 
Electric Field 


It is seen from above that in order to get the solutions of 
the electric intensity and electric displacement, we can 
turn to solving the modal potential functions. So, the 
modal Equation (39) should have the corresponding mo- 
dal boundary condition. 

An effective boundary case is: Electric displacement 
functions of two side of interface should be equal 


"=" so 


or 


z 
7 


D? i=1,2,3 (47) 


i i 


Rewriting Equation (46) in the modal form, we have 


ae a 


or 
AMM yO = AC ye% i=1,2,3 (49) 
and 
AGH = AVY" i=1,2,3 (50) 
7. Application 


In this section, we discuss the laws of static electric field 
only in anisotropic dielectrics. 


7.1. Isotropic Crystal 


The matrix of dielectric permittivity of isotropic dielec- 
trics is following 


é&, 0 0 
é€=|0 «, 0 6D 
0 0 «&, 


The eigen-values and eigen-vectors are respectively 
shown as below 


A= diag|&,,&,,€| (52) 
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100 
®=|0 1 0 (53) 
001 


We can see from the above equations that there is only 
one eigen-space in isotropic crystal, which is a triple- 
degenerate one, and the space structure is following 


W =W [404] (54) 


The basic vector of one dimension in a triple-degenerate 
subspace is 


8 


t= {LLD (55) 


The eigen electric displacement operator of isotropic 
crystal are 


z.i 
vi = 5 (2 +4; +85) (56) 


a NB 
{1} =F 10.22.83} (57) 
Therefore, the static electric field equation in isotropic 


crystal can be written as below 
(02 +02 +07) x, =0 (58) 


Thus, the electric strength and electric displacement of 
isotropic crystal become 


1 
{E} = {4} Vizi = (8) +03 +83) l} z 
1 (59) 
1 
{D} = {Bj AVinn =e, (a? +05 + ð; ) l% 
0) 
The classical electric potential function of order 1 is 
yı =—(8, +0, +83) Z; (61) 


So, the electric intensity and electric displacement of 
isotropic crystal can also be expressed by the classical 
electric potential function of order | as follows 


ô, 
(E}=-42, fy (62) 
ô, 
or 
E =-Vy, (63) 
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ô, 
{D} STE ô, VY, (64) 
ô, 
or 
D = -6 Vy, (65) 


It is seen that Equations (62)-(65) are the same as the 


classical results, in which the electric potential is a scalar. 


But from the following analysis, we will see that only for 
isotropy we have same results as classical theory. 


7.2. Uniaxial Crystal 


The matrix of dielectric permittivity of uniaxial dielec- 
trics is following 


ér 9 0 
é=|0 & 0 (66) 
0 0 &, 


The eigen-values and eigen-vectors are respectively 
shown as below 


A =diag[|£1; £183] (67) 
10 0 

®=|0 1 0 (68) 
001 


We can see from the above equations that there are 
two eigen-spaces in uniaxial crystal, one of which is a 
double-degenerate space, and the space structure is fol- 
lowing 


W =W [gp] m [g] (69) 
The basic vectors in two subspaces are following 


¢ = ei 1,0% (70) 


g, ={0,0,1)" (71) 


The eigen electric strength qualities of uniaxial crystal 
are respectively shown as below 


E; =) -E=E, (72) 

ġ Ei =E-$E, (73) 

Multiplying Equation (57) with %, using ¢) -¢, =0 
and g'-g =1 (i=1,2), we get 


= (2-8) (E-@°E})= JE? +E} (74) 


The eigen electric displacement operators of uniaxial 
crystal are respectively shown as below 


* 


1 


E 
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Vi =5(8u+dm) (75) 
V} = ô% (76) 

a Ne 
{a} = “> {01,82,0} (77) 
{a} = {0,0,8,} (78) 


Therefore, the static electric field equation in uniaxial 
crystal can be written as below 


(22 +22) zi =0 (79) 


ox =0 (80) 


It is seen from Equations (79) and (80) that there are 
two static electric fields in uniaxial crystal. Thus, the 
electric intensity and electric displacement become 


{E} = {a} Viti +{¢,} eee 


1 0 
=;1 (gO) + 07 033% (81) 
0 1 
(A, +055 vA 
= (ô + 3z rA 
One 
{D} z {Oh AVin +{b}AVin 
£E (ô tea) (82) 
= £ (ôa +3a) mh 
£33033% 


The classical electric potential function of order 1 is 
1 
=-—(0,+0 
WY, 2 ( 1 2 xm (83) 
yi =% (84) 


So, the electric intensity and electric displacement of 
uniaxial crystal can also be expressed by the classical 
electric potential function of order 1 as follows 


OW, 

{E}=-4d,y, (85) 
OW, 
E0 N 

{D} = -5 818 V (86) 
E330.YW> 


It is seen that the electric intensity and electric dis- 
placement of uniaxial crystal are quite different from 
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those in isotropic crystal, and there exist two kinds of 
modal electric potential functions in uniaxial crystal, so 
they become vectorial ones, this is also different from the 
classical results of static electric field. 


7.3. Biaxial Crystal 


The matrix of dielectric permittivity of biaxial dielectrics 
is following 


é, 0 0 
e=|0 e, 0 (87) 
0 0 £; 


The eigen-values and eigen-vectors are respectively 
shown as below 


A = diag | £1; E2, E3] (88) 
100 

ð=0 1 0 (89) 
001 


We can see from the above equations that there are 
three eigen-spaces in biaxial crystal, and the space struc- 
ture is following 


W =w. [OM [6] OM” [g] (0) 


The eigen-qualities and eigen electric displacement 
operators of biaxial crystal are respectively shown as 
below 


Ei =¢, -E=E, (91) 
E,=$,-E=E, (92) 
E; =$;-E=E, (93) 


Vi =4, {0} ={8,,0,0}" (94) 
Vi=0y  {0,}={0,0,,0)' (95) 


V;=0, {0;}={0,0,0,}" Wo 


Therefore, the static electric field equation in biaxial 
crystal can be written as below 


aixi =0 (97) 
oix = 0 (98) 
atx; =0 (99) 


It is seen from Equations (97)-(99) that there are three 
static electric fields in biaxial crystal. Thus, the electric 
intensity and electric displacement become 
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{E} = {01} Viz +9} V22 +13} V32 


1 0 0 
=40}V 7, +41 bV +0% (100) 
0 0 1 
Onn 
= 40% 
033%3 
EON 
{D} = Enaka (101) 
£33033- X3 


The classical electric potential function of order 1 is 


Yı =-0,% (102) 
WY, =-05% (103) 
V, = -03% (104) 


So, the electric intensity and electric displacement of 
uniaxial crystal can also be expressed by the classical 
electric potential function of order 1 as follows 


OV 
{E} = -40V 
OW; 


(105) 


RORA 
{D} == Ey 0 Wo 
E330_W, 


(106) 


It is seen that the electric intensity and electric dis- 
placement of biaxial crystal are quite different from 
those in isotropic crystal, and there exist three kinds of 
modal electric potential functions in uniaxial crystal, so 
they become vectorial ones, this is also different from the 
classical results of static electric field. 


7.4. Monoclinic Crystal 


The matrix of dielectric permittivity of monoclinic di- 
electrics is following 

E E2 0 

E=| E2 En 0 (107) 


0 O &, 


The eigen-values and eigen-vectors are respectively 
shown as below 


A = diag | 4,4, £] (108) 
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T 
$, = biz aa 1 o 
1 z 3 6 z oi 
-6,) +E 12 
(A, i) 12 (109) 
$ za i h -£n 0 
2 5 Rie) ax 
(A, - 4) +E 12 
g, =[0, 0,1] 
where 
2 
(4: tEn) 1 2 
Ay = 2 =f Me Ey) +€) A, =&; (110) 


We can see from the above equations that there are 
also three eigen-spaces in monoclinic crystal, and the 
space structure is following 


W =W [a] om,” [a ]Om [g] 


The eigen-qualities and eigen electric displacement 
operators of monoclinic crystal are respectively shown as 
below 


(111) 


1 
E =ġ -E= 7 [(4-81)E, +£E, | 
A, £1) +E 
(112) 
* 1 
E, =): = [EE +(4-61) | 
(A, —&,) +E 
(113) 
E; =¢,-E=E, (114) 
V? =a20,, +b23,, +2a,b,ô,ð, i 
{0} = {a,ð, +,0,,4,0, +b,2,,0}" 
V? = 420, + b28, +2a,b,0,0, aie 
{a5 | = {a,0, +b,0,,,0, +b,0,,0}" 
V;=0,, {03} ={0,0,0,}" (117) 
where 
b, = fa a, =) j=1,2 
(4 =en) + E15 E2 


Therefore, the static electric field equation in mono- 
clinic crystal can be written as below 


(a70,, +70, +2a,)0,0,) 0 (118) 
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(a30,, +B20,. +2a,b,0,0,) 73 =0 (119) 


aty, =0 (120) 


It is seen from Equations (118)-(120) that there exist 
also three static electric fields in monoclinic crystal, 
which is a little different from the results in biaxial crys- 
tal because of the distortion of static electric fields. Thus, 
the electric intensity and electric displacement become 


{E} ={¢, bv 141+ {o}V SA +g} Vi = 
aĝ +60, 
að, +b,0, (a7d,, + Bf ô» +2a,b,0,0,) x, 


4,0, a b,0, 
ı +5,0, (a30,, +b), + 2a,b,6,0, ye 


(121) 
0 


003343 
1 


+ 


{D} z {Ot AVin TAZAS +{Ø} 4V% T 
a,0,+ 50, 
4,0, +b, 7A (adi + bi Oy + 24,b,0,0, eA 
0 
a,0, +b,0, 

+} 4,0, +b,0, tA, (a30, + Ay, + 2a,b,0,0,) x, 
0 
0 


+3 0? 4033% 
1 


(122) 


The classical electric potential function of order 1 is 


v= -[a (4,0, +b0,)+b, (a0, +b0, )] a (123) 


W= -[a, (a, +b,0,)+b, (4,0, +b,0, )] Xa (124) 


W, = -0; %3 (125) 


So, the electric intensity and electric displacement of 
monoclinic crystal can also be expressed by the classical 
electric potential function of order 1 as follows 


aĝ, + b,0, 4,0, + b,0, 0 

{E} =—4 4,0, +b,0, py, =) 4,0, + b,0, pW —40 PW, 
0 0 ô; 

(126) 
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a0, +50, 
{D} = -4 440, + 5,0, pW 
0 
a,0,+b,0, 0 
—A, 44,0, +b,0 (Wy — E3359 FY; 
0 0; 


(127) 


It is seen that the electric intensity and electric dis- 
placement of biaxial crystal are aiso quite different from 
those in isotropic crystal. 


8. Conclusions 


In this paper, we construct the standard spaces under the 
physical presentation by solving the eigen-value problem 
of the matrixes of dielectric permittivity and magnetic 
permeability, in which we get the eigen dielectric per- 
mittivity and eigen magnetic permeability, and the cor- 
responding eigen vectors. The former are coordinate- 
independent and the latter are coordinate-dependent. Be- 
cause the eigen vectors show the principal directions of 
electromagnetic media, they can be used as standard 
spaces. Based on the spaces, we get the modal equations 
of static electromagnetic fields by converting the classi- 
cal Maxwell’s vector equation to the eigen Maxwell’s 
scalar equation, each of which shows the existence of an 
static electromagnetic field. For example, there is only 
one kind of static electromagnetic field in isotropic crys- 
tal, which is identical with the classical result; there are 
two kinds of static electromagnetic fields in uniaxial 
crystal; three kinds of static electromagnetic fields in 
biaxial crystal and three kinds of distorted static electro- 
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magnetic fields in monoclinic crystal. All of these new 
theoretical results need to be proved by experiments in 
the future. 
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